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Abstract 

Resonance trapping appears in open many-particle quantum systems at high level density when 
the coupling to the continuum of decay channels reaches a critical strength. Here a reorganization 
of the system takes place and a separation of different time scales appears. We investigate it under 
the influence of additional weakly coupled channels as well as by taking into account the real part 
of the coupling term between system and continuum. We observe a saturation of the mean width 
of the trapped states. Also the decay rates saturate as a function of the coupling strength. The 
mechanism of the saturation is studied in detail. In any case, the critical region of reorganization 
is enlarged. When the transmission coefficients for the different channels are different, the width 
distribution is broadened as compared to a Xk distribution where K is the number of channels. 
Resonance trapping takes place before the broad state overlaps regions beyond the extension of the 
spectrum of the closed system. 



1 Introduction 

At low excitation energy, the states of a many-particle quantum system are usually well isolated from 
one another. Their coupling via the continuum of decay channels can therefore be neglected. The 
reaction cross section consists of a sum of resonances with Breit-Wigner shapes the positions, widths 
and partial widths of which are well defined. 



At higher excitation energy both the level density and the decay widths of the states in most sys- 
tems become so large that the resonance states overlap. The cross section is an interference picture. 
Neither the positions, the widths nor the partial widths of the resonance states can unambiguously be 
determined from an analysis of the cross section, see e.g. [Q. Meaningful values can be obtained from 
time measurements which provide the average lifetime of the resonance states at high level density. 
Using the channeling method, the mean lifetime of fine structure resonances under isobaric analogue 
resonances in medium nuclei was found experimentally to be much longer than expected on the basis 
of the statistical theory of nuclear reactions . The modification of the Fano profiles of autoionizing 
states due to the coherent coupling with each other has been studied experimentally by means of the 
two states 3p^ and 3p3d in magnesium atoms Q 



Theoretically, the transition from low to high level density is studied recently in different papers 
for different systems, see |]|] and |Q to [26|. As a function ofT/D (where T is the average value of the 
widths of all resonance states, being a measure for the coupling strength of the system to the contin- 
uum, and D is the mean level distance) all results show the same characteristic features: beyond a 
critical value ofV/D, separated time scales exist when the number of channels is smaller than the num- 
ber of resonance states. The very existence of different time scales at high level density corresponds to 
the basic assumption of the unified theory of nuclear reactions formulated phenomenologically about 
40 years ago by Feshbach for many resonances coupled to a small number of common decay chan- 
nels. 



The main differences in the results of the different theoretical approaches consist in the behavior 
of the long-lived resonance states as a function of further increased T/D. While the widths saturate 
in the more realistic models for many-body systems e.g. [14|, they decrease in the random matrix 
models. A saturation of the decay rates is observed also in the "bottle-neck" picture of transition state 
theory which relates the saturation value to the number of independent decay channels rather than 
to the widths of individual resonance states of the system |I2S 



As a reply to |29|, the mechanism of 
the saturation is shown by the same authors to be associated with a broadening of the distribution of 
resonance widths [^]. A broadening of the width distribution at high level density is found also, e.g., 

to be caused by resonance trapping which is the basic 



in p|, pi and shown in |, W^, M, pT 



process of the redistribution of the system. It leads ultimately to the formation of different time scales. 



In this paper, we study in detail the widths of the long-lived states as a function of the coupling 
strength to the continuum since this is a controversial point of discussion. Most calculations are per- 
formed beyond the standard random matrix approach. 

In sect. 2, the phenomenon of resonance trapping is described in the framework of the random 
matrix theory. The widths of the trapped resonance states do not saturate with increasing coupling 
strength to the continuum. In sect. 3, the spectroscopic values of an open many-particle quantum 
system are given. Additonal terms in the effective Hamiltonian appear which may prevent the widths 
of the trapped resonance states to decrease with increasing coupling to the continuum. In the following 
sections, we discuss the influence of additonal terms in the effective Hamiltonian on the widths of the 
long-lived trapped states. Additional weakly coupled channels are shown to cause a saturation of the 
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mean width with further increasing couphng strength (sect. 4). The sharp distinction between short- 
and long-hved states is removed by the real part of the couphng term to the continuum (sect. 5). In 
both cases, further avoided crossings appear and the bi-orthogonahty of the eigenfunctions of H^^ 
continues to be essential. That means, the critical region where the redistribution takes place is en- 
larged. This can be seen also in the width distribution of the trapped states, investigated in sect. 6. 
Some conclusions on the broadening of the width distribution, the separation of time scales and the 
saturation of decay rates at high level density are drawn in the last section. 



2 Random matrix theory and resonance trapping 
2.1 Basic equations 

In order to describe resonance phenomena at high level density in a small energy interval far from 
thresholds the random matrix theory has been developed. The effective Hamiltonian is (see e.g. 

M i I) 



^RMM = H^-iW = H^- '-VV^ . (1) 

Here V is an energy independent random matrix consisting of K random vectors of dimension with 
matrix elements V^. The value = Sr=i(^)^ gives the coupling strength of the system to the 
channel c. Each vector Vc has Gaussian distributed elements with mean value and variance /N . 
Hf, is chosen from the Gaussian Orthogonal Ensemble (GOE) and the mean level density follows a 
semi-circle law. In our case the length of the spectrum is L = 2 units and therefore the mean level 
distance in the middle of the spectrum is given hjd = tt/{2N). 



According to 33 1, one can obtain the mean width (F;) of the long-lived resonance states by 
considering the diagonal elements of the S-matrix, averaged over a sufficiently large energy interval. 
It holds 



K 



X{\{S,c)\=e--^^^^'' . (2) 

c=l 

The mean width (T^) of the states in the middle of the spectrum is 

(r^) = -^EMi-rc), (3) 



27r . 

c=l 



where the Tc are the transmission coefficients, defined as follows, 

r, = l-K5,,)|2. (4) 



These coefficients are calculated in |34] using a resummation method in the power series expansion of 
the transfer matrix. The result is 

_ 4xe _ ttI^ , . 

- (1 + ^ J2 ' - • l^j 

The mean width T of all resonance states can be obtained directly from (||), 
f = llm {tr <mm} = ^EK^ = • 
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The parameters Xc measure the couphng strengths of the system to the different decay channels c. For 
increasing total coupling they subsequently pass over the critical value one. Here the corresponding 
transmission coefficient reaches its maximal value Tc = 1. The logarithm in (^) diverges and Eq. (|3|) 
looses its validity. Eq. (|5|) is symmetric with respect to an exchange of rcc by 1/xc. As a consequence 
the transmission coefficients do not distinguish between weak and strong total coupling strength. We 
have Tc < 1 for all Xc 7^ 1. 

In the following sections we will see, that at each of the points Tc = 1 a single resonance state 
separates in width from the remaining ones. For Xc > 1 this state becomes much broader than the 
other ones and can no longer be described statistically together with the long-lived resonance states. 



The width of the short-lived state is |11, 25 1 



2Nd 
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Xc > 1 (7) 



while the widths of the long-lived states decrease as l/x^. 



2.2 Resonance trapping at strong coupling 

One interesting question is how the system described by the hamiltonian H^^^y^ behaves as the cou- 
pling to the continuum increases. To investigate this question, we replace V by aV in the term W of 
the effective Hamiltonian, Eq. (||), and vary a. 

For small a we can treat W as a, small perturbation on Hi). The rank of //f, is N. Hj^j^jj^ thus 
describes states with energies determined by the eigenvalues of Hf,. The widths are proportional to 
a^. This holds well as long as the resonances are non-overlapping. 

For large a we can treat Hi, as a small perturbation on W. As W has rank K it directly follows 
that only K states have large widths, see e.g. 0]. The widths of the remaining N — K states are 
small. That means, two different time-scales exist. 

In this subsection, we first illustrate the phenomenon of resonance trapping for A^ resonance states 
coupled to K = 1 open decay channel. We define 



N 



Here the sum runs over all but the j broadest states. Dq = L/{N — 1) is the mean level spacing of the 
eigenvalues of Hi, and L ist the length of the spectrum. For N — j ^ 1 the difference between k^~^ 
and is 2Tj/{DQ{N — j)) (states ordered according to decreasing width). Further we define 

tot _ ,j=o _ A V — - 2— (()) 



K 

M 

R=l 



K*"* is a measure of the total coupling strength of the system to the continuum (see Eq. (^)). In the 
random matrix theory (with A^ large) we have d = 7r/(2N) and L = 2 and thus Dq = Ad/ir which 



gives K*"* = Xc in the one-channel case (compare subsection |2.lD . 

Within the random matrix theory, we calculate versus for A^ = 300, K = 1 and j = 0, 1, 2 
with varying coupling to the continuum (fig. l.a). The curves shown are averages over 20 calculations. 
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Note the logarithmic scales. Until ~ 1 the average width of all states increases with increasing 
coupling to the continuum. At «:*°* ~ 1 two globally separated time-scales are formed. The broad 



state should be identified with a doorway state [13| and a tight transition state |35|, respectively. At 



the separation point = 1, the transmission coefficient is Tc = 1. For still further increasing 
the broadest state is getting a still larger width but the average width of the remaining ones decreases. 

The thick line in fig. l.a. shows 2(J'i)/Dq obtained from Eqs. (|3|) and (^). For k*"' < f, (F;) 
is the mean width of all resonance states, (F;) w F, while for At*°* > 1, (F;) is the mean width of 
the A — 1 trapped states, (F^) ~ Fjv^i. Therefore, the transmission coefficients Tc related to the (F/) 
by Eq. (H) do not give us information on the coupling strength of the system to the continuum 



(compare subsection 2.1) 



Resonance trapping of many states takes place whenever the local level density compared to the 
local mean width is sufficiently large, see e.g. |2l[ ]. 

To illustrate the local properties of resonance trapping we show in fig. l.b the widths F/j versus 
the energies for the states of fig. l.a. Only a small part of the spectrum is shown. The calculations 
are performed for 0.1 < < 10 in steps of log«;*°* = 0.04. The points for = 0.1, 1 and 10 
are marked with triangles, stars and squares, respectively. As a function of a the complex eigenvalue 
er = Er — of each resonance state follows a certain "trajectory". For small coupling to the 
continuum, the widths of the states increase with increasing continuum coupling. This process takes 
place for every resonance state up to that value of the continuum coupling at which the state starts to 
overlap one of the resonance states in its neighborhood. The crossing of resonance states is avoided in 
the complex plane: The states attract each other in energy and their widths bifurcate, i.e. the width 
of one of the states continues to grow with further increasing coupling strength while that of the other 
one decreases. The state finally being the broadest one goes through a number of "collisions" (avoided 
resonance crossing between two states looks like a "collision" in this representation) before it dominates 
the complete spectrum. It is formed in the middle of the spectrum where the level density is the largest. 

The stars on the trajectories (fig. l.b) show that trapping of resonance states is a local process 
which takes place before different timescales are formed globally. The widths of all the trapped states 
decrease with further increasing coupling strength. 

This behaviour does not change when the level density has a band or shell structure. We simulate 
such a situation by choosing a H\^, Eq. (|l]), with several more or less separated regions of high level 
density. Between them and at the borders of the spectrum, the level density is smaller. In the present 
calculation we have K = A and A = 300. The coupling of each state to each channel is randomly 
chosen, i.e. the coupling strengths of the four channels are equal to one another. We have three regions 
with Gaussian shaped level density. Two of the regions are close to each other and the third one is 
lying well separated from the other two. The energies E^ and widths of the resonance states are 
shown in fig. 2 for varying total coupling strength to the continuum. 

We see the following result. In each group of states there are formed some broad states at high 
level density. These states do not overlap regions outside the interval studied. With further increasing 
coupling strength, these broad states attract one another in energy and their widths bifurcate. Finally, 
there are only four broad states according to the four channels. In the separation point, their widths 
are smaller than the energy region covered by the long-lived states. The question whether there are 
states lying outside the extension of the spectrum is of interest only when studying the behaviour of 
the short-lived states at a further increased strength of the coupling to the continuum. The resonance 
trapping in every group is left unchanged by the effects at the borders of the spectrum. 
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We conclude this section by stating the following. In the random matrix theory, the widths of the 
trapped resonance states decrease as a function of increasing coupling strength k*"*. 



3 Continuum shell model (CSM) for an open many-body quantum 
system 

In many-body quantum systems, the widths of the long-lived resonance states do not decrease but 
saturate as a function of the coupling strength 3C]. In order to investigate the mechanism of 



saturation we have, according to the results shown in subsection 2.2, to go beyond the random matrix 



theory. In the following, we sketch a model which allows to describe open many-body quantum systems. 

The time independent Schrodinger equation 

{H-E)\'^) = (10) 

is solved in a Hilbert space consisting not only of the discrete many-particle states of a closed system 
but also of the continuum of decay channels. The potential is assumed to be a spherical one. The 
Hamilton operator of the system is 

H = H^ + V. (11) 

is the unperturbed Hamilton operator describing particles in a finite depth potential and V is the 
operator of the two-body residual interaction between the particles. For details see [^. The relation 
between V , Eq. (|l|), and the two-body operator V is considered in |36|. 



In order to find the solution of (jl^) we first solve the shell model problem 

[E^ - Fqq) I = (12) 

in the Q subspace of N discrete states and the coupled channel equations 

(i^(+)_i/pp)|^^(+)) = (13) 

with the proper boundary conditions in the P subspace of K coupled channels. The projection 
operators are Q = J2r=i I ^T)i^T I and P = Ef=i I dE \ \ . Here i/qq = QHQ and H^^ = 

PHP. Further, we solve the coupled channel equations with source term 

(i?(+)-//pp)i4+)) = ypQi<^fr>, (14) 

which connects the two subspaces. Using P -|- Q = 1 we then express |^') by means of the solutions 
of (^^, (13) and (0). Care must be taken in order to avoid double-counting from ( |l^ and (13), 



i.e. appearance of any resonances in (M). For this purpose a cut-off technique for single particle 
resonances is used in when solving (ll^). 

In P subspace, the propagator is G^^^ = P{E + ie — Hpp)~^P . The propagator Q{E — Hqq}~^Q 
in Q subspace contains the effective Hamilton operator in this subspace, 

H'^(:^iE) = Hqq + Vqp4^^Vpq . (15) 
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H^q{E) is non-hermitean and energy dependent. It has energy dependent complex eigenvalues 
Er = Efi — ^Tr and eigenfunctions | ^r) describing the quasi-bound states embedded in the con- 
tinuum (QBSEC) 1^. The | ^r) form a bi-orthogonal set at each energy E with the orthogonality 
relation {^rI^'r) = 5rr' (see Further, it holds {^r\^r) > 1 , {^r\^'r) eC {R^R') . 



Diagonalizing -f^gq we get the solution of (10) as 



N . 

1 4^^) =1 ct^) + E I ^^R^)-^{^n I y I d^^) • (16) 

R=i ^-^R 

Here the 

0^^) = {1 + G^p^^V)\^r) (17) 

are the wavefunctions of the resonance states R. It should be stressed here once more that V is the 
operator of the two-body residual interaction and that the ^r are many-particle wavefunctions. 

The relation ( |T7| ) between the wavefunctions ^^r^ of the resonance states and the eigenfunctions 
^R of H-^Q is analo gous to the Lippman-Schwinger equation 

lei) = (1 + 4^V) ixi) (18) 

between channel wavefunctions Xe ^^'^ coupled channel wavefunctions Therefore {Qr\V\x^^^^) = 
. We define the amplitude of the partial width by 

7i?,e = -^mVlxt^) = -±={^R\V\ct^) . (19) 



By inserting the expression (16) for ^' into the S-matrix and using ( [19| ) we get 

R E - Er-\- -^L R 

The first two terms in ( pC| ) describe the direct reaction part of the process. The last term describes the 
resonance part, i.e. excitation and de-excitation of the resonance states R. The resonance part S^f^ 
of the S— matrix has the standard Breit-Wigner form. Note however that ^rc is complex and energy 
dependent. Also Er and Tr are energy dependent functions. S'^^f^{E) contains the contributions of all 
the resonance states at the energy E of the system. These contributions may be very different from 
those at other energies (see [p^]). 



The poles of the S-matrix give the energies Er and widths Tr of the resonance states. They are 
defined by the fixed-point equations Er = Er{E = Er) and Tr = Tr{E = Er). The number of 
resonance states is exactly equal to the number of discrete states obtained from (|l^) if double 
counting in P and Q is avoided. 

If the energy dependence of Er{E) and Tr {E) is weak in the interval studied, it holds Er ~ Er{E') 
and Tr ~ Tr{E') where E' is an energy somewhere in the middle of the interval. Only under such 
conditions, Er and Tr can be considered as resonance parameters. Otherwise, the spectroscopic stud- 
ies performed at the energy E' are meaningful only at this energy since only here the orthogonality 
relations between the right and left wavefunctions are fulfilled. The present calculations are performed 
at the energy E = E' oi the system. 
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According to [^], the decay rates can be defined by = ln((/>(t)|(/)(t)) where (j){t) is the 

wave function of the system. Using the ansatz \(j){t)) = J2r o^Rit) \^r) and solving the time-dependent 
Schrodinger equation with the effective Hamiltonian -f^qq we get 

\m = E aij(0)e-i(^«-5f «) * 1$^) . (21) 
R 

The aji{0) define the wavefunction <p(0) of the system at the time t = 0. Neglecting the oscillations 
caused by the bi-orthogonality of the function system we get 

in analogy to = Tr/Ti for isolated resonances. The are given by = |aij,(0)p 

g-^Rt/^ ($^1$;^) . They decrease exponentially with the rate Tr/Ti. Therefore the for the 

short-lived states are neglible in the long-time scale and the sum in (^2[) runs only over the trapped 
states. If the width distribution of the N — K trapped states is narrow, it holds f ~ T]^^k- The 
weighted width T*^"^) is generally time dependent. 



H^SoiE) = Hqq + V{E) - iW{E) . (23) 



Let us write Eq. (|15|) more explicitly 

jeff 



Here 

wrr, (s) = ^ E mmR' (E) (24) 

c 

where V^{E) = {(l)j^\V\xE'^^) are real numbers describing the coupling of the shell model states R to 
the channels c at the energy E of the system and c runs over all open channels. Further, 

oo 

Vrr,{E) =Y.Pj VRiE')^^V^,{E')dE' . (25) 

P denotes the principal value of the integral and Cc is the threshold energy for the channel c. As a 
rule, the €c are different for different channels. 

Characteristic of the Hamiltonian of an open many-body quantum system is therefore (i) the differ- 
ent channels are coupled with different strength and (ii) the coupling via the continuum contains not 
only an imaginary part but also a real part. In the following sections we study the widths of trapped 
resonance states as a function of increasing coupling strength by considering these properties of 
the Hamiltonian. 



4 Influence of additional weakly coupled channels 

In a many-body quantum system, each channel has a certain coupling strength which may be quite 
different from the coupling strength of other channels. Sources for the different coupling strengths of 
the channels are, above all, the structure of the states of the residual system, the different angular 
momenta and the different threshold energies €c, see sect. 3. 
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We study the influence of different coupling strengths of the channels onto the resonance trap- 
ping in a schematical manner. The main emphasis lies on the question whether the mean width of 
the trapped states at strong coupling to the continuum saturates or approaches zero. The study is 
performed in the random matrix model (RMM) with the Hamiltonian (||) and N = 300, K = 4, but 
different average coupling strength to the channels. The ratios among the coupling strengths are 
V^/V^^i = 1, 0.1, 0.01 and 0.001 (fig. 3.a) and V^/V^^i = 1, 0.01, 0.0032 and 0.001 (fig. 3.b). We 
show for j = 0...5 averaged over 20 calculations. 

The separation points are defined by Tc = 1 for any channel, i.e. Xc = 1 for a certain c (see Eqs. 
(1) and (D). These points are given by = J2d=i^d/^c- This gives = 1.11, 11.1, 111, 1110 
in fig. 3.a and = 1.01, 101, 321, 1010 in fig. 3.b. 

In fig. 3. a, there is at ~ 1.2 a separation point, where one broad state separates from the other 
ones and the mean width of the remaining resonance states starts to decrease. Here, for j > 
decreases slightly, but soon increases again under the influence of the next channel. At ~ 12 a sec- 
ond broad state separates from the other ones. Similar situations occur at k*"* ~ 120 and ~ 1200. 
Finally there are four broad states corresponding to the four channels. For still larger values of k*"*, 
the values , j > 3 decrease. Note that the values , j > 4 are almost constant in the region 
1 < < 1500. Thus, different coupling strengths to the channels are a source for saturation of the 
average width of the long-lived N — K states as a function of . 

In fig. 3.b, we have one strongly coupled channel and a group of weakly coupled ones. At k*"* ~ 1.1 
a picture similar to ordinary trapping with one channel can be seen. One state separates and the 
widths of the remaining ones decrease. For larger coupling to the continuum however, the widths of 
the trapped states start to increase under the influence of the weakly coupled channels. A new critical 
region occurs between « 70 and 2000 where three states corresponding to the three new channels 
separate. In this region, the widths of the trapped states saturate and thereafter their widths decrease. 

Further, these two examples show another interesting result. At the points where a broad state 
separates from the remaining ones, we have Tc ~ 1 for a certain channel c, see the thick line {2{Ti)/Do 
versus When all Tc < 1, we have good agreement between the (F/) calculated from Eqs. (|3|) and 

(Q) and the calculated mean width Tn-j- Here, j is determined by the number of broad states and 
j = only for < 1 (compare subsection |2.2| ) . 

A smah part of the eigenvalue picture {Er and Tr), fig. 3.c, (V^/V^^i = 1, 0.01, 0.0032 and 0.001 
as in fig. 3.b, but only one calculation) shows width increase, energy shift, width decrease and once 
again increase, shift and decrease. The renewed increase of the widths of the states trapped by the 
first channel is caused by the fact that the new channels become active only at strong coupling to the 
continuum. 

In fig. 3.d, we show Tr/2 and {xc — l/xc)/2 , c = 1,...,4 , Eq. (|^), versus for the same 
calculation as in fig. 3.c. The agreement between (xc — l/xc) and the four different widths F^ of the 
short-lived states illustrates nicely that the separation of every broad mode takes place from the group 
of long-lived resonance states. This means, the separation process is more or less independent of the 
broad modes separated at smaller values of 



5 The role of the real part of the coupling term 
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5.1 Two resonances coupled to one common channel 



The influence of a real part onto the phenomenon of resonance trapping can be seen by means of the 
fohowing simple model for two resonance states coupled to one common decay channel, 

h^''=(l ^)-2^ae-^'( ^"I'^f" ^ 1 ; "90° < ^ < 90° . (26) 

\ — 1 / \ cos (p sm if sm 9? / 

Here the relative coupling strength of the two states to the continuum may be varied by means of the 
angle ip {ip ^ 0'^,90'', ...). The angle 9 determines the ratio between the real and imaginary part of 
the coupling term. In ll^ the case 9 = has been studied. 



The eigenvalues of h*^^ are 

= -iae"*^ ± ^1 - 2iae-^^ cos{2lp) - a'^e-'^^^ . (27) 

At a = we have two states lying at the energies —1 and 1. Their widths increase with increasing a 
up to some critical value Ocrit ^ 1- The eigenvectors of h'^^ are 

$i = ' { \ . (28) 



Here, the normalization of the wavefunctions is made according to (^*±*|*^±) = 1 and i?i>± = 1 — 
2ia e~*^ cos^{ip) — £± .In our case of only two resonance states, the values |$±P are the same for 
both states. 



The distance in the complex plane between the two eigenvalues is 

2 



|e+ — e_ 



1 — 2iae *^ cos(2<^) — a^e 



2S 



Using S, the denominator in Eq. (p8|) can be rewritten as follows. 



a2 e-2^e sin2(2(^) + 0^ = J2S (S ± (1 - ie"^^ cos(299))) . 



(29) 



(30) 



The factor {S it (1 — *^ cos(2(^))) is never zero for finite S, and thus |<I>p oo if and only if the 
distance in the complex plane between the two eigenvalues goes to zero. 



Let us define k = (Fi + T2)/L where L = 2 is the distance between the two resonance states. For 
the hamiltonian (^^ it holds k = 2 a cos 9. 

In the following the two cases f = 22.5^ and = 45'^ are studied in detail. In the first case both 
states coupled with equal strength to the decay channel while in the second case the state having the 
initial energy +1 is stronger coupled to the channel than the other one. 

In fig. 4 the complex eigenvalues (4. a) and the I'I'p versus k (4.c) are shown for ip = 45^ and 9 = 0^, 
W^, 30° and 45°. The different points in fig. 4. a corespond to different values of k, 0.2 < n < 20. As k 
grows, each eigenvalue follows a certain trajectory. For = 0° the Eji and F/j of the two states meet 
in one point in the complex plane at «; = Kcrit = 2. Here |^*(Kcrit)P oo. Beyond this separation 
point, one state continues to increase in width whereas the width of the other one decreases {resonance 
trapping). For 9 > the state with initial energy —1 becomes the broader one and is shifted towards 
negative energies. In contrast to this, for all < (not shown in the figures) the state with inital 
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energy +1 becomes, in this symmetrical case 93 = 45 , the broader one and is shifted towards large 
positive energies. For 9 ^ the minimum distance in the complex plane between the two states 
remains different from zero and |$| remains finite for all k with its maximum value at /^crit • 

As 

9 90*^, h^^ becomes hermitian. In this limit, |<I>p — > 1 and Kcrit 0. 

The case = 22.5*^ is shown in figs. 4.b (complex eigenvalues) and 4.d versus k) for some 
values of 9 ranging from —45" to +6O". For ^ > the state with initial energy = +1, being the 
broader one at small k, gets an extra shift towards small energies whereas the shift is towards large 
energies for < 0. Kcrit is a function of only \9\ and has the same values as in the symmetrical case. It 
has its largest value for 9 = . At 9 = Afp the complex eigenvalues of the two states meet in one point 
and |<I>(Kcrit)P diverges. ks9 ^ ±90° h'^^ becomes hermitean and |<I)(k™*)P 1. For -90° < 61 < 45° 
the state with initial energy = — 1 is the one becoming trapped. For 45° < 9 < 90°, however, the 
state with initial energy Eji = +1 becomes trapped even though that state is the broader one at small 

K. 

These examples show that the details of the resonance trapping change when allowing for extra 
energy shifts by introducing the angle 9 in (p6|). They are basic for an understanding of the results of 
the following section. 

5.2 N resonance states coupled to K common channels 

We study the influence of the real part of the coupling term onto the mean width of the trapped states 
in the framework of the CSM. We do this by comparing the results of calculations with and without 
V, Eq. (|2^), taken into account. 

The rank of "P is, generally, larger than K. V increases with a by approximatively the same factor 
as W. Thus the hermitean part Hqq + V of is never a small perturbation on W and the rank of 
i/qq is larger than K also in the strong coupling limit (compare sect. |2.2D . Nevertheless, resonance 
trapping occurs also when taking V into account. 

In fig. 5 we present versus for some calculations in the CSM. As in the RMM, we use 
Dq = L/{N — 1) where L is the length of the spectrum of Hqq. We study 190 resonance states in ^^O 
with J'^ = 1~ and vary a. In the first case (figs. 5. a, b) we study the reaction ^^Oi/2-{n',ny^Oi/2- 
giving K = 2 channels with s and d waves, respectively. The second case (figs. 5.c, d) is {Pi 
where x = p, n and Y = 3/2-5 ^^Oi/2-, 3/2- which gives K = 10 channels. Figs. 5. a, c and 

5.b, d are without and with V, respectively. 

The figures 5. a, c with "P = are similar to figs. 3. a, b. In fig. 5. a we have 2 channels with 
different coupling strength. First we see the increase of the mean width of all the states for small 
coupling to the continuum up to the separation point corresponding to the strong channel. Thereafter 
we have a region of saturation of the mean width and the separation point for the second broad state. 
In fig. 5.C we see first the separation of 4 strongly coupled states and the saturation of the widths up 
to the coupling strength at which 10 states corresponding to the 10 channels are separated from the 
remaining ones. Finally the mean width of all the trapped states decreases. 

In the calculations with the principal value integral V taken into account, we additionally see 
another effect (figs. 5.b, d). There are, at strong coupling to the continuum, not only K broad states, 
but more states separate from the remaining ones and get large widths. In fig. 5.d we have at = 60 
about 15 states that are broad but the sharp distinction between the lifetimes of the broad and the 
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trapped states is washed out. 

Similar results are obtained by studying other resonance states with other quantum numbers. 

In the eigenvalue picture, fig. 6. a, we show and Tji for the same resonance states as in fig. 
5.d by varying in the interval 0.008 < < 60. Note that the steps in are approximatively 
equividistant. We see extra shifts in energy caused by the principal value integral V. Such a shift is, 
generally, in the order of magnitude of the width of the state, i.e. it is large for states getting large 
widths with increasing coupling strength to the continuum. That means, the broad states leave the en- 
ergy region where the trapped states are lying. These shifts are similar to those in figs. 4. a, b for 6^0. 

Fig. 6. a illustrates also the behavior of trapped states under the influence of increasing coupling 
to the continuum. First their widths increase, then the states get trapped, i.e. their widths start to 
decrease and they get a small energy shift. For even stronger coupling to the continuum the widths 
can start to increase again with a renewed energy shift and so on. In distinction to figs, l.b and 3.c, 
these shifts have mainly two origins: the energy attraction accompanying the bifurcation of the widths 
and the influence of V. 

The corresponding as a function of are shown in fig. 6.b. The |^_rP remain larger than 1 

also at large values of k*"*. This indicates further avoided resonance crossings which cause, ultimately, 
the saturation of the widths of the long-lived resonance states. 

Summarizing, the V enlarges the critical region of reorganization where the local process of reso- 
nance trapping takes place. The mean width of the long-lived states saturates but there is no longer 
a sharp distinction between the broad states and the long-lived ones. 



6 Resonance trapping and broadening of the width distribution 



In the RMM with equally strongly coupled channels it is shown in e.g. [26, ^ that the width distri- 
bution in the critical region of reorganization is broader than in non-critical regions. In non-critical 
regions, i.e. as long as all transmission coefficients are small, the width distribution follows a Xk la-w 
with the number of degrees of freedom corresponding to the number K of open channels Q. In this 
section we study the width distribution in the RMM with varying coupling strengths to the channels 
and in the CSM. 

In the RMM with varying coupling strengths of the channels we have tried to fit the width dis- 
trubution in non-critical regions to a Xr distribution. The error in the fit is large even for small 
transmission coefficients. We conclude that the system must be coupled to all the channels with com- 
parable strengths for the widths to be Xk distributed. 

In the following we therefore study the broadening of the width distribution by calculating the 
normalized variance of the widths. 



1 

\ ^— E - 1)' yn = r«/f (31) 

\ R=l 



where the sum runs over all but the j broadest states. The theoretical value obtained in the RMM 
with K equally strongly coupled channels for the (trapped) states far from the critical region of reor- 
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ganization is ai=^ = y^J/K = af^^. 

The results of calculations in the RMM with varying coupling strengths to the channels are shown 
in fig. 7. In fig 7. a a-jj is shown for K = 10 channels with V^/V^^i distributed on [1 ... 0.1] with equal 
distances at the logarithmic scale for j = 0...13. In figs. 7.b and c we show for K = A channels 
and j = 1...6 with V^/V^^i = 1, 0.1, 0.01, 0.001 and V^/V^^i = 1, 0.01, 0.0032, 0.001, respectively 
(compare figs. 3. a, b). Note that cr^ before the separation of a broad state should be compared with 
the value cr^"''^ after the separation. 

In the figures, the values a^^"^ = ^IjK are shown with a dashed line. In the case of different 
coupling strengths of the K channels, the width distribution is broader than in the case of K channels 
with comparable coupling strengths. 

The separation of every broad state is accompanied by a broadening of the width distribution, see 
figs. 7. a to c. The distribution between the separation points in fig. 7.b and between the separation 
points of the first and second broad state in fig. 7.c. is however narrower than the distributions for 
very small and very large k*"*. 

To describe the distribution of the long-lived states in regions where all Tc < 1 we make the 
following ansatz, compare Eq. (^: 

r« = -^Ei^(i-^c)5L- (32) 

c 

Here ^ for a certain c is a Gaussian distributed vector with mean zero and variance 1 and is the 
transmission coefficient defined in Eq. (^). 

The normalized variance of is 



, V ("^^^^ - v/2Ee (ln(l-r.))^ 

= FT = -Eeln(l-r.) ^-'y ■ 



(Note that = \ and = 3.) ay = a^^^ for the same number K of open decay channels holds only 
when all are equal. The values calculated from ( |33|) are shown with a thick line in figs. 7. a to c. 
The variances are well described by (^) in regions of the coupling strength with all Tc < 1 . 

The values as a function of k calculated in the CSM are similar to those obtained in the RMM. 
One example is shown in fig. 7.d for the case J'^ = 2^ with K = 12 open channels. 

We conclude that in many-body quantum systems, the distribution of the widths is broader than 
in the RMM with equal coupling strength to the channels (also at small transmission coefficients). 
This result explains the broadening of the width distribution described in the literature by introducing 
an effective number K^^ of channels being smaller than K. 

For the case shown in fig. 7. a, we have calculated also k^, Eq. (p2|), as a function of time t with 
12 different coupling strengths k*"* of the system to the channels (fig. 8). h = 1. In table Q the values 
^tot g^^^ Ti^-K for the different curves are given. All the curves a to f, being below the critical region, 
differ considerably from one another. Inside the critical region, the curves for different (g to i) 
are, however, similar to one another. 
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As a result, weakly coupled channels cause a saturation of the mean value Tn-j as a function of 
K^°^. Since the width distribution does not change much, also saturates in the long-time scale, 
i.e. remains almost unchanged by varying in the critical region. The saturation is related to 
the enlarged width distribution. 



7 Conclusion 

In this paper we studied the positions and widths of resonance states in a many-body quantum system 
with N resonance states as a function of increasing coupling strength k*"' to the continuum which 
consists K <^ N decay channels. In a critical region of the coupling strength, the system reorga- 
nizes itself under the influence of the decay channels. The local process is the avoided crossing of two 
resonance states which takes place whenever the distance in energy between the states is comparable 
to the sum of their widths. It is accompanied by an essential bi-orthogonality of the eigenfunctions 
of H'^^. As a result, one state continues to increase in width whereas the other one decreases with 
further increasing coupling strength (resonance trapping). For a system with many states, this leads 
to a broadening of the width distribution. 

With further increased coupling strength, the broadening of the width distribution goes over into 
a separation of time scales if K <^ N. The smaller K the better expressed is this separation. Weakly 
coupled channels as well as the hermitean part of the coupling term VqpGp'^VpQ have the tendency 
of washing out the differences between the lifetimes of the group of long-lived trapped states, on the 
one hand, and the group of short-lived states, on the other hand. 

We studied in detail the widths of the N — j long-lived trapped states under different conditions. 
We introduced additional channels to which the system is weakly coupled and we took into account 
the hermitean part of VqpGp^^VpQ in Hqq, Eq. (15). Under the influence of these additional terms 



in the Hamiltonian, the trapped states can increase their widths and change their positions in energy. 
Thus, trapped resonance states may again come close to each other if the coupling to the continuum 
is stronger. The resonance crossing is avoided and accompanied by an essential bi-orthogonality of 
the eigenfunctions of H^^ in the same manner as at smaller coupling strength. As a result of all 
these processes, the average width of the states saturates as a function of the coupling strength to the 
continuum when an appropriate number of states with the largest widths is excluded from the mean 
value. This number is equal to the number of open decay channels as long as the real part of the 
coupling term VqpGp^^VpQ is small compared to its imaginary part. 

The width distribution of the long-lived states is related to the transmission coefficients. If all 
transmission coefficients are equal and smaller than one, the widths are Xk distributed where K is 
the number of channels. In many-body quantum systems with different coupling strengths to the 
different decay channels, the distribution is broader than a Xk distribution also for small transmission 
coefficients. 

The decay rates are related to the mean decay widths of the long-lived states. This means, also the 
decay rates saturate in quantum systems at high level density as a function of the coupling strength 

„tot 

Summarizing the results we state that resonance trapping is a realistic process occuring in many- 
particle quantum systems at high level density. It leads to a saturation of both the decay rates and the 
average decay width of the long-lived states as well as to a broadening of the width distribution and 
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- if the number of open decay channels is not too large - to a separation of time scales. That means, 
the decay rates and the decay widths of the long-lived resonance states show the same behavior not 
only at low level density but also at high level density. The saturation is caused in both cases by the 
finite number K < N of channels into which the N resonance states can decay. The application of the 
standard random matrix approach to the details of the trapping process in an ensemble of resonance 
states is limited. 

We would like to state once more that the time de-excitation of resonance states at high level 
density should be directly measured. The results could make a proof of the phenomenon of resonance 
trapping possible. 
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Figure 1 



for j = 0,1,2 and 2{Ti)/Dq, calculated from Eq. (thick line), versus (l-^)- Eigenvalue 
picture {^Tr and Er) calculated for different k*"' (l.b). The calculation shown is performed within 
the RMM for K = 1 and N = 300. In (l.a) the curves shown are averages over 20 calculations. Only 
a part of the spectrum is shown in (l.b). Note the logarithmic scales. The points for = 0.1, 1 and 
10 are marked in (l.b) with triangles, stars and squares, respectively. 



Figure 2 

Eigenvalue picture (^F/j and Er) for H\y consisting of a sum of three Gaussian shapes. = 300, K = 4 
and all channels are coupled with the same strength. 



Figure 3 

for j = 0...6 versus in the RMM for N = 300 and K = 4 with the ratios of the coupling 
strengths Yc/^c=i = 1> 0.1, 0.01, 0.001 (3.a) and V^/V^^i = 1, 0.01, 0.0032, 0.001 (S.b). Eigen- 
value picture (3.c) and ^Fr (dots) together with the function ^{xc — 1/xc) , Eq. (|^, (full line) versus 
(3.d) for V^/V^^i = 1, 0.01, 0.0032 and 0.001. In (3.a) and (3.b), 2(r,)/L>o versus is shown 
with a thick line. In (3.c), the points for = 0.1, 1, 10, 100 and 1000 are marked with triangles, 
stars, diamonds, large dots and squares, respectively. 



Figure 4 

Eigenvalue picture (^E/j and Er) calculated for different k (4. a, b) and |^>p versus k (4.c, d) for 
two states coupled to one channel. It is (/? = 45° (4.a, c) and tp = 22.5° (4.b, d). The values are shown 
for some different 6 in the same plots (see the text for details). Note the logarithmic scales. 



Figure 5 

versus in the CSM for 190 resonance states in ^^O with J'^ = 1~ , K = 2 (5. a, b) and K = 10 
(5.C, d) channels. In (5. a, c) "P = while V is taken into account in (5.b, d). It is j = 0...4 (5. a, b), 
j = 0...14 (5.c) and j = 0...25 (5.d). The calculations are performed at £^ = 29 MeV. 



Figure 6 

Eigenvalue picture {^Fr and Er) for different coupling strengths 0.008 < < 60 (6. a) and I^r]"^ 
versus k}°^ (^-b) in the CSM. N = 190 resonance states in ^^O with J'^ = 1~ and = 10. V is 
taken into account (compare fig. 5.d). 
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Figure 7 



(7^, ay and a^^^ versus in the RMM for three selected values of V^/V^^^ (for details see the 
text), is shown for j = 0...14 (7. a) and j = 0...6 (7.b, c). The thick lines are ay obtained from 
Eq. @) and the dashed hues are a^^'^^. (7.d): a^ for j = 1..20 versus in the CSM for J"" = 2", 
K = 12 and V is taken into account, j = 0...20 and the dashed line is a^^^ for K = 12. 



Figure 8 

k„ versus time t for if = 10 and V^/V^^i = 1...0.1 (compare fig 7.a). The 12 curves are calculated 
with between 0.1 and 50. and Fat-x for the curves are given in table ffl. 



Table 1: k*"* and T^-k for the different curves in fig. 8 



Curve Tn-k 



a 


0.10 


0.094 


b 


0.18 


0.17 


c 


0.31 


0.29 


d 


0.54 


0.51 


e 


0.96 


0.89 


f 


1.7 


1.5 


g 


3.0 


2.6 


h 


5.2 


3.9 


i 


9.2 


4.7 


j 


16. 


4.9 


k 


28. 


4.5 


1 


50. 


3.3 
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